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Abstract 

Ammunition  deterioration  during  storage  has  considerable  eco- 
nomic consequences.  A  reliable  prediction  model  for  the  ammunition 
deterioration  rate  is  necessary  for  long-term  procurement  and  main- 
tenance planning.  A  random  effect  growth  curve  analysis  is  employed 
to  formulate  a  prediction  model  for  ammunition  deterioration  rates  in 
terms  of  concurrent  characteristics  such  as  depot  condition  and  vendor 
information.  The  resultant  prediction  model  can  be  used  to  determine 
the  appropriate  time  for  reorder  or  renovation  of  ammunition  before 
the  quality  reaches  unacceptable  levels.  A  two-stage  analysis  is  used 
to  estimate  parameters  involved  in  the  prediction  model.  Necessary 
estimation  methods  are  discussed.  An  example  is  given  to  illustrate 
the  implementation  procedure  of  the  prediction  model  suggested  in 
this  paper. 
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1.  INTRODUCTION 

Each  year,  the  U.S.  military  services  purchase  millions  of  dollars  worth 
of  major  caliber  ammunition  from  several  vendors  to  prepare  for  possible 
war.  Ammunition  lots  purchased  are  stored  in  depots  until  they  are  used.  In 
view  of  the  fact  that  most  ammunition  is  produced  long  before  its  ultimate 
consumption,  it  is  important  that  the  material  be  adequately  stored  and 
maintained  to  remain  usable. 

However,  according  to  the  GAO  report  (The  Comptroller  General  (1982)), 
the  military  services  are  experiencing  significant  problems  in  maintenance 
of  depot  and  renovation  of  ammunition.  These  problems  cause  ammunition 
to  deteriorate  faster  than  necessary  maintenance  can  be  performed  and  a 
serious  backlog  of  ammunition  in  need  of  maintenance.  In  an  attempt  to 
prevent  such  undesirable  phenomena,  in  this  paper,  a  quantitative  model 
that  can  relate  a  depot  maintenance  plan  to  the  quality  control  of  ammu- 
nition is  considered. 

The  objective  of  this  paper  is  to  outline  a  prediction  model  for  ammuni- 
tion deterioration  in  terms  of  depot  characteristics  and  vendor  information. 
The  goal  is  to  provide  estimation  methods  by  which  one  can  determine  the 
appropriate  time  for  reorder  or  renovation  of  ammunition  before  the  qual- 
ity reaches  unacceptable  levels.  The  suggested  model  can  also  be  used  for 
vendor  control  and  examination  of  depot  maintenance  policy. 

In  order  to  formulate  such  a  prediction  model  for  ammunition  dete- 
rioration rates,  a  random  effect  growth  curve  analysis  is  employed.  The 
analysis  of  a  random  effect  growth  curve  model  has  been  extensively  used 
in  biomedical  research  (Rao  (1965),  Stiratelli  et  al.  (1984),  and  Vonesh 
and  Carter  (1987)).    An  application  of  growth  curve  analysis  to  physical 


engineering  is  given  in  Sohn  and  Mazumdar  (1991).  In  general,  the  goal  of 
growth  curve  analysis  is  to  estimate  random  growth  rates  of  a  time  varying 
principal  attribute  in  terms  of  concurrent  individual  characteristics  which 
do  not  change  over  the  experimental  period.  The  analysis  of  growth  curves 
is  typically  carried  out  by  using  a  two-stage  random  effect  model.  The 
first  stage  consists  of  the  within-individual  regression  model  in  which  the 
within-individual  growth  rates  for  the  serial  observations  are  estimated  for 
all  individuals.  In  the  second  stage,  estimated  growth  rates  are  related  to 
a  set  of  covariates  representing  concurrent  individual  characteristics. 

In  this  paper,  a  growing  pattern  of  defective  proportion  of  an  attribute 
of  ammunition  lot  over  time  (deterioration  rate)  is  viewed  in  analogy  to 
a  growth  rate  of  a  principal  attribute  of  an  individual  experimental  unit 
while  vendor  information  and  depot  characteristic.-  associated  with  each 
lot  can  be  considered  as  concurrent  individual  characteristics.  In  order  to 
predict  deterioration  rates  of  a  certain  type  of  ammunition  in  terms  of  depot 
characteristics  and  vendor  information,  a  random  effect  logistic  regression 
model  (Stiratelli  et  al.  (1984))  is  employed  assuming  that  the  pattern  of 
deterioration  follows  a  logistic  function.  For  estimation  of  parameters  in 
a  random  effect  logistic  regression  model,  a  two-stage  estimation  based  on 
the  method  used  in  Korn  and  Whittemore  (1979)  is  employed.  An  example 
is  given  to  illustrate  the  implementation  procedure  of  the  random  effect 
logistic  model  suggested. 


2.  RELATED  STUDIES 

The  most  important  deterioration  processes  occurring  to  ammunition 
during  storage  are  corrosion  of  the  metallic  components  and  chemical  and 
physical  changes  within  the  internal  structure  of  the  explosives,  which  may 
result  in  reduced  stability  and  performance  (Eckman  (1974),  Eriksen  and 
Str0ms0e  (1978),  and  Wilson  and  Bird  (1986)).  The  factors  which  have  the 
greatest  influence  on  these  processes  are  temperature  and  humidity  which 
are  directly  related  to  the  location  and  condition  of  the  depot.  Even  under 
similar  conditions,  different  types  of  ammunition-  may  degrade  with  differ- 
ent deterioration  rates.  In  some  instances,  different  vendors  provide  the 
same  types  of  ammunition  that,  nonetheless,  have  varying  deterioration 
rates.  Most  of  the  studies  in  the  literature  do  not  reveal  the  quantita- 
tive relationship  between  deterioration  rates  and  depot  condition,  types  of 
ammunition  and  vendor  information. 

For  instance,  Eriksen  and  Stromsoe  (1978)  conducted  an  extensive  ex- 
periment covering  a  9  year  period  to  examine  deterioration  of  six  types  of 
artillery  ammunition  (105  mm  HE,  90  mm  HE,  57  mm  HE,  40  mm  HE, 
hand  grenades  Mk2,  and  antitank  mines  M6A2)  at  three  aboveground  and 
three  underground  depot  sites.  Various  attributes  were  observed:  vacu- 
um stability  and  water  content  of  propellant,  impact  sensitivity  of  primers, 
burning  rate  of  propellant,  muzzle  velocity,  and  the  number  of  misfires  and 
duds.  Inspection  frequencies  varied  from  two  to  four  times  covering  com- 
binations of  starting  year  (1968)  and  the  subsequent  1st,  5th,  7th  and  9th 
year  after  storage,  depending  upon  the  attributes  inspected.  Based  on  these 
repeated  measurements,  the  contingency  table  and  the  analysis  of  variance 
were  used  to  test  if  the  levels  of  ammunition  deterioration  were  significantly 


different  in  terms  of  the  conditions  of  depot  and  the  types  of  ammunition. 
The  results  of  their  study  are  as  follows:  Overall,  aboveground  sites 
surrounded  by  inland  climate  performed  better  than  underground  depots 
that  were  old  and  technically  in  poor  shape;  the  visual  inspection  revealed  a 
considerable  difference  in  corrosion  attack  between  the  ammunition  types; 
water  content  both  in  the  primer  and  the  propellant  charges  varied  signif- 
icantly depending  on  the  depot  in  which  lots  of  ammunition  were  stored; 
most  of  the  chemical  analyses  were  of  little  interest  and  should  be  postponed 
to  the  end  of  trials.  Their  approach,  however,  did  not  provide  a  compre- 
hensive model  which  related  the  level  of  deterioration  of  certain  type  of 
ammunition  to  depot  condition.  If  the  level  of  deterioration  over  time  can 
be  expressed  as  a  function  of  depot  conditions  or/and  vendor  information, 
reorder  time  of  a  certain  type  of  ammunition  which  reaches  certain  quality 
can  be  predicted  in  terms  of  given  information  at  the  depot.  An  appli- 
cation of  growth  curve  analysis  which  facilitates  such  a  prognostication  is 
illustrated  in  the  following  section. 


3.  A  GROWTH  CURVE  ANALYSIS  FOR  AMMUNITION 
DETERIORATION  RATE 

Consider  a  certain  type  of  ammunition  (say,  105mm  HE).  In  this  section, 
the  assumptions  used  in  the  growth  curve  analysis  to  predict  the  deteriora- 
tion rate  of  105mm  HE  ammunition  are  summarized.  Similar  assumptions 
can  be  made  for  other  types  of  ammunition. 

Basic  unit  of  purchase  is  a  lot  and  as  a  result  of  acceptance  sampling, 
the  qualities  of  incoming  ammunition  lots  would  be  homogeneous  regardless 
of  different  sources  of  manufacturers.  While  several  hundreds  of  such  lots 
are  stored  in  various  locations,  the  quality  of  ammunition  deteriorates  with 
the  passage  of  time.  Deterioration  rates  may  vary  depending  upon  both  the 
condition  of  depot  in  which  ammunition  lots  are  stored  and  manufacturers 
who  provide  the  ammunition  lots. 

In  order  to  inspect  the  pattern  of  deterioration  of  various  lots  over  the 
time,  it  is  assumed  that  a  set  of  initially  selected  lots  in  various  depot- 
s  is  marked  and  these  lots  are  repeatedly  inspected  on  a  sampling  basis 
without  rectification  over  a  certain  time  period.  Lot  size  is  assumed  to 
be  sufficiently  large  that  it  would  not  be  depleted  at  the  end  of  repeated 
sampling  inspections.  Modes  of  inspection  can  cover  visual  inspection,  lab- 
oratory examination  and  functional  trials  (Eriksen  and  Stromsoe  (1978)). 
The  proportion  of  samples  which  do  not  satisfy  a  predetermined  level  of  an 
attribute  or  the  combinations  of  attributes  can  be  used  as  a  proxy  for  the 
accumulated  level  of  deterioration  at  the  time  of  inspection.  For  instance, 
among  many  other  attributes  of  functional  trials,  the  relative  number  of 
misfires  and  duds  to  the  sample  size  observed  at  each  inspection  cycle  can 
be  considered  as  a  measure  for  the  level  of  deterioration  of  each  lot  at  the 
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time  of  inspection. 

For  each  lot,  the  proportion  of  defective  rounds  is  assumed  to  increase 
as  time  goes  on  following  a  logistic  curve.  The  logistic  model  is  well  es- 
tablished as  a  basis  for  analyzing  such  a  phenomenon  when  the  response 
variable  follows  binomial  distribution.  The  logistic  model  formed  to  de- 
scribe the  deteriorating  pattern  for  each  lot  is  called  a  within-individual 
model.  Deterioration  rate  of  a  lot  of  ammunition,  which  represents  a  dete- 
riorating pattern  over  time  of  a  lot  often  differs  from  that  of  another  lot. 
The  model  in  which  varying  deterioration  rates  between  lots  are  compared 
in  relation  to  depot  characteristics  and  vendor  information  associated  with 
each  lot  is  called  the  between-individual  regression  model.  Depot  condition- 
s  and  vendor  information  associated  with  each  lot  of  certain  ammunition 
are  assumed  to  remain  constant  during  the  experimental  period  while  they 
would  vary  over  different  lots.  Necessary  notations  for  a  two-stage  random 
effect  logistic  model  for  ammunition  deterioration  rate  is  as  follows: 

NOTATIONS 

i:  index  for  individual  lot,  i  =  1, ..,  N 

N:  total  number  of  ammunition  lots  employed  in  the  experiment 

j:  index  for  inspection,  j  =  1, ..,  nt 

n,-:  the  number  of  inspections  to  be  done  on  a  lot  i  over  experimental 
period 

tij-.  jth  inspection  time  of  a  lot  i 

Yij\  sample  size  of  lot  i  at  the  jth  inspection 

Xijk'-  a  dichotomous  variable  indicating  if  each  item  k  taken  from  the 
sample  at  the  jth.  inspection  of  lot  i  is  defective  or  not 

y 

Xij\  the  number  of  defectives  found  in  a  sample  size  of  >£;  Xl3=Y^k=\  xijh 


Pq  =  Xij/  Yij'.  proportion  defective  at  the  jth  inspection  of  a  lot  i 

exp(j30) / (l+exp(j30)):  initial  proportion  defective  of  ammunition  lot 

/3t:  the  deterioration  rate  of  lot  i  (log  odds  ratio  over  the  unit  time 
interval  for  lot  i) 

(2,1, 2i2,  ■■yZim):  covariates  indicating  concurrent  characteristics  of  lot  i 
(e.g.,  dummy  variables  indicating  different  combinations  of  depot  condition 
and  vendor  information  of  ammunition  lot  i) 

71, ..,  7m:  regression  coefficients  for  (2,1, ..,  zim) 

m:  total  number  of  covariates  in  the  between-individual  model 

7:  a  pxl  vector  of  7's,  (71, ..,  jm)' 

Z:  an  Nxm  matrix  of  covariates  (2,1, ..,  2im) 

e,:  random  error  oilog(3t  ~  N(0,  cr2). 

p*:    estimated  expected  proportion  defective  of  lot  i  at  time   t  given 

%:  estimated  expected  time  that  quality  of  lot  i  reaches  a  predetermined 
level  of  proportion  defective  p  given  (z&,  ..,2tm) 

/?o,  Pi'-  maximum  likelihood  estimators  of /30  and  (3X 

log(3:  an  iVxl  matrix  of /o^/3,'s,  (log(3\,..,log(3Ny 

Si',  estimation  error  associated  with  log(3t,  6t  ~  N(0,  if ) 

7q,  t?:  variance  of  j3q  and  that  of /3, 

7^,  tf:  maximum  likelihood  estimators  of  7^  and  rt2 

W:  an  NxN  diagonal  matrix  with  diagonal  elements,  l/(o"2  +  rf  //?2)'s 

7:  maximum  likelihood  estimator  for  7 

Mathematical  formulations  regarding  the  random  effect  logistic  model 
are  as  follows: 

In  the  within-individual  model,  proportion  defective  of  lot  i  at  the  j th 

8 


inspection,  p^  is  assumed  to  follow  a  logistic  function  of  time  t^: 
For  i  =  1, ..,  N,  and  j  —  1, ..,  ?i, 

p#  =  exp(pQ+Pitij)/{\  +  expiPQ+PiUj))  (1) 

or 

log(pi,/(l-pij))=p0  +  (3itii  (2) 

In  the  between-individual  model,  deterioration  rate  of  ammunition  lot 
i,  Pi,  is  represented  as  a  function  of  covariates  (e.g.,  depot  condition,  and 
vendor  information)  as  well  as  a  random  error  (e,): 

Pi  =  fillZii  +  ■■  +  ImZim  +  e.)  (3) 

where  e,  follows  independent  N(0,  a2).  The  proportion  defective  of  ammu- 
nition lot  would  not  decrease  as  time  goes  on.  The  following  function  is 
used  to  ensure  positive  deterioration  rates,  i.e.,  positive  /3,'s: 

Pi  =  exp(jizn  +  ..  +  ymZim  +  e»)  (4) 

or 

log(J3i)  =  71-2,1  +  ••  +  JmZim  +  C,  (5) 

After  replacing  /?,■  in  (2)  with  that  in  (4),  one  could  estimate  regression 
parameters  (/?0,  7i,--,7m5  and  cr2),  by  maximizing  the  marginal  likelihood 
of  the  data  based  on  the  combined  model.  However,  exact  solutions  for  the 
normal  equations  based  on  the  combined  model  are  analytically  and  com- 
putationally infeasible  (Stiratelli  et  al.  (1984)).  Indicating  these  problems, 
Stiratelli  et  al.  (1984)  discussed  alternative  ways  of  estimating  parameters 
(EM  algorithm  and  a  two  stage  estimation).  When  P,  in  (2)  is  replaced  with 


that  in  (4),  random  error  e.'s  appear  as  coefficients  of  the  regression  model. 
In  EM  algorithm  they  are  treated  as  missing  observations  and  are  esti- 
mated by  cycling  back  and  forth  between  the  E-step  and  the  M-step  until 
convergence  is  reached.  When  a  two-stage  estimation  method  is  used,  the 
estimation  procedure  for  the  within-individual  model  is  separated  from  that 
for  the  between-individual  model.  For  instance,  in  the  within-individual  re- 
gression model,  deterioration  rates  of  a  set  of  individual  ammunition  lots 
are  estimated  for  all  experimental  lots.  In  the  between-individual  regression 
model,  the  estimated  deterioration  rates  are  used  as  values  of  a  dependent 
variable  and  they  are  regressed  on  depot  characteristics  and  vendor  infor- 
mation associated  with  each  individual  lot. 

In  general,  a  two-stage  estimation  method  is  computationally  less  inten- 
sive and  analytically  easier  than  the  EM  algorithm  (Stiratelli  et  al.  (1984)). 
In  the  following  section,  a  two-stage  estimation  is  described  in  order  to  for- 
mulate a  prediction  model  for  ammunition  deterioration. 
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4.  A  TWO-STAGE  ESTIMATION 

First,  in  order  to  estimate  deterioration  rates  in  the  within- individual 
model  (1),  the  following  likelihood  function  ofz^  conditional  on/30,  Pi,..,Pn 
is  formulated: 


N     n,     X,  N     n,     /  y 

V 

x 


l = n n  n pw  - ^)l_i" =nn?  k**1  - **>*~x'    (6) 

v=lj  =  lk=l  i=\j  =  l\^V/ 


where  p^  is  as  in  equation  (1). 

Maximum  likelihood  (ML)  estimates,  /?o,  P\,..,/3n  are  obtained  by  solv- 
ing a  set  of  normal  equations  of  log  likelihood  function  of  (6)  with  respect 
to  /?0,  /3i,..,(3N.  Estimated  parameters  (/3,-'s)  are  in  turn  assumed  to  be  ap- 
proximately normal  with  means  equal  to  the  true  individual  parameters 
(/3,'s)  and  variances  (if 's).  Diagonal  elements  of  the  inverse  of  the  nega- 
tive information  matrix  evaluated  at  (flo,  /3i,..,/3at)  provide  the  estimated 
variances  of  (/?o,  ft,..,^),  which  are  (tq  ,  ff, ..,  fjy). 

Secondly,  resultant  /?,-  replaces  unobservable  /3,  in  the  between- individual 
model.  This  replacement,  however,  brings  the  estimation  error  <5,  to  the 
equation  (5): 

\og0i)  =  JiZa,  +  ..  +  ImZim  +  Ei  +  St.  (7) 

where  6t  is  statistically  independent  of  e,.  It  is  assumed  that  <5,  follows  in- 
dependent Normal  distribution  with  mean  0  and  variance  if  /j32  (var  (<$,-)  = 
var(/o^/3,  |  /?,•)  ~  var  ((/?,-/?  ,•)//?,•)  =  t2/(32).  The  var(<5l)  is  approximated  by 
•ff //??  and  /og/3,  is  marginally  distributed  as  independent  N(712ll  +  ,  ■•,+7m2tm, 
a1  +  r1 /fl2).  Standard  normal  theory  can  be  used  to  estimate  the  between- 
individual  model  parameters  tj's  and  a2  which  maximize  the  following  log 
likelihood  function  at  the  estimated  r2  : 
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N  N  m 

LL  =  -Y,log(a2  +  if)  -  ^(logfc  -  £  7^)7(^2  +  «?//§?)•        (8) 

»=i  1=1  1=1 

The  resultant  MLE  of  7  is 

y=(Z'WZ)-\Z'W(log/3))  (9) 

and  the  estimated  variance  of  7  is  (Z'WZ)~  where  an  NxN  diagonal 
matrix,  W,  consists  of  diagonal  elements,  1/(<t2  +  if  //?2)'s.  In  order  to 
be  able  to  estimate  7,  the  number  of  lots  (TV)  employed  in  the  experiment 
must  be  lager  than  the  number  of  covariates  (m)  to  be  used  in  the  between- 
individual  model.  In  addition,  as  N  increases,  accuracy  related  to  the 
inference  concerning  7  would  increase. 

When  these  7/s  replace  7/s  in  (7),  logfti  can  be  predicted  in  terms  of 

log(/?J  =  7l2»l  +  ••  +  Trn^trn-  (10) 

Use  of  (10)  to  quality  control  of  ammunition,  vendor  control  and  depot 
maintenance  plan  is  illustrated  in  the  next  section  along  with  a  numerical 
example. 
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5.  IMPLIMENTATION 

Unfortunately,  actual  data  related  to  type,  quantity,  and  depot  condi- 
tion of  ammunition  were  omitted  from  the  GAO  report  (The  Comptroller 
General  (1982))  to  keep  it  unclassified.  In  the  Norwegian  study  cited  in  the 
earlier  part  of  this  paper,  repeated  measurements  on  the  deterioration  of 
ammunition  was  not  exposed  for  all  the  experimental  lots.  In  fact,  only  the 
summary  statistics  for  deterioration  of  several  groups  of  ammunition  were 
available,  respectively.  In  order  to  illustrate  implementation  procedures  for 
the  methods  suggested  in  this  paper,  a  numerical  example  is  given  based  on 
the  simulated  data  generated  partly  based  on  the  parameters  used  in  the 
Norwegian  study  (Eriksen  and  Stromsoe  (1974))  and  the  guideline  of  the 
U.S.  Army  Ammunition  Surveillance  Procedure  (Supply  Bulletin  SB  742-1 
(1988)). 

Table  1  contains  information  regarding  the  series  of  the  number  of  defec- 
tive items  (Yjj)  found  in  20  ammunition  lots  (TV).  These  lots  are  classified 
into  four  groups  depending  upon  the  vendor  who  provided  ammunition  and 
depot  in  which  ammunition  is  stored:  group  1  consists  of  5  lots  (i  —  1, ..,  5) 
which  are  provided  by  vendor  A  and  stored  in  depot  1;  group  2  consists  of 
5  lots  (i  =  5, ..,  10)  which  are  provided  by  vendor  B  and  stored  in  depot  1; 
group  3  consists  of  5  lots  (i  =  11,  ..,15)  which  are  provided  by  vendor  A 
and  stored  in  depot  2;  and  group  4  consists  of  5  lots  (i  =  16, ..,  20)  which 
are  provided  by  vendor  B  and  stored  in  depot  2.  All  the  lots  are  inspected 
annually  based  on  sample  size  of  20  (Y,j).  Total  number  of  inspections  of 
each  lot  (n.)  varies  over  4  different  groups:  n,  =  7  for  i  —  1, ..,  5;  n,  =  8  for 
i  =  6, ..,  10;  n,  =  9  for  2  =  11, ..,  15;  and  m  =  9  for  i  =  16, ..,  20. 

Our  goals  are  first  to  find  the  relationship  between  deterioration  rates 
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and  concurrent  characteristics  of  ammunition  lot  (depot  condition  and  ven- 
dor information);  secondly  to  determine  if  any  group  of  ammunition  causes 
significantly  different  deterioration  rates  from  the  others;  and  finally  to  pre- 
dict the  time  for  renovation  of  ammunition  in  terms  of  characteristics  of 
depot  and  vendor. 

In  order  to  achieve  the  first  goal,  a  two-stage  estimation  is  employed. 
First  of  all,  the  ML  estimates  for  parameters  in  the  within-individual  model 
(/?o,  /3i,.., /?2o  and  f02,  r^,..,^)  are  obtained  based  on  (6).  For  this  step, 
PROC  LOGISTIC  of  a  statistical  package  SAS  (SAS  (1989))  is  used  by 
applying  the  values  of  XtJ,  }£,  t^  in  Table  2  to  model  (2).  Resultant  /3(s 
and  rf's  are  summarized  in  Table  2  while  /30  turns  out  to  be  -4.9034  which 
indicates  that  the  estimated  initial  proportion  defective  of  ammunition, 
(exp(j30)/(l  +  exp0o))),  is  about  0.74%  regardless  of  sources  of  vendor. 
Sample  patterns  of  the  actual  deterioration  ( Actual:  A"^/  1^  verses  t^)  are 
given  in  Figure  1  along  with  those  of  the  estimated  deterioration  (Fitted: 
(exp0o  +0itg)/(l  +  exp0o+l3lttj))  verses  fy). 

The  between-individual  model  is  formed  to  relate  deterioration  rates  to 
the  ammunition  characteristics.  log/3t  is  used  as  dependent  variable  and  the 
four  dummy  variables  (z&, ..,2*4)  representing  4  combinations  of  depot  and 
vendor  characteristics  are  used  as  covariates  without  an  intercept:  2,1  =  1 
for  ammunition  lots  which  belong  to  group  1,  2,-1  =  0  otherwise;  Z&  =  1  for 
ammunition  lots  in  group  2,  zq  =  0  otherwise;  z^  =  1  for  ammunition  lots 
in  group  3,  z&  =  0  otherwise;  2,4  =  1  for  ammunition  lots  in  group  4,  z^  =  0 
otherwise.  Given  such  (2,1,  ..,z&)  as  well  as  (/?-,  and  rf's),  an  optimization 
package  GAMS  (GAMS,  1988)  is  applied  to  (8)  in  order  to  obtain  71, ..,  74 
and  a1.  Estimates  are  summarized  in  Table  2  along  with  standard  errors  of 
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7i, ..,  74  obtained  from  the  square  root  of  the  diagonal  elements  of  (Z'  WZ)   l 

(U      0      0      0\ 


where  a  20x4  covariate  matrix  Z  = 


;  15  =(1,1,1,1,1)'  is 


0      15      0      0 

0      0      15      0 

V0      0      0      U) 
a  5x1  unit  vector;  and  Wis  a  20x20  diagonal  matrix  consisting  of  diagonal 

elements  1/(<j2  +  f,-2//?2).  Values  of  r2  and  /32/s  are  given  in  Table  1  while 

(72=0.012. 

Estimated  7  (I  =  1,  ..,4)  is  the  weighted  average  of  logfi'fi  which  belong 

to  the  group  /.    A  large  7  indicates  fast  deterioration.    For  instance,  the 

estimated  average  log  deterioration  rate  of  ammunition  lots  in  group  2  (- 

0.8336)  is  lager  than  that  of  ammunition  lots  in  group  1  (-0.9022).  To  test 

if  any  group  characteristic  causes  significantly  faster  deterioration  rate  than 

another,  a  (1  —  a)  100%  confidence  interval  for  Z'7  can  be  used: 


l'j±  t(a/2,N  -p)y/l'(Z'WZ)-H 


(11) 


where  /  is  a  4x1  vector  which  contrasts  a  pair  of  two  group  characteristics; 
t  is  a  student  t  distribution  with  degrees  of  freedom  N  —  p. 

In  order  to  test  if  the  difference  in  two  expected  log  deterioration  rates 
between  group  1  and  group  2  is  significant,  one  can  use  /'  =  (1,  —1,0,0)  to 
obtain  a  95%  confidence  interval  for  ji  —  72: 

(1,-1,0,  0)(-0.9022,  -0.8336,  -0.7184,  -0.6567)' 


±<(0.025, 20  -  4)v/(l,-l,0,0)(Z'fyZ)-1(l, -1,0,0)' 

/  0.00708  0  0 


where  (Z'WZ)'1  = 


\ 


0 
0 
0 


0.00473  0 

0  0.00278 


0 


0 


0 
0 
0 
0.00226/ 
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The  resultant  interval  (-0.0686  ±  2.12(0.10870))  contains  0  which  im- 
plies that  at  a  =  5%  vendor  effect  is  not  significant  on  the  deterioration  of 
ammunition  stored  in  depot  1.  Similar  comparisons  can  be  made  by  mod- 
ifying I's.  Complete  analyses  of  95%  confidence  intervals  for  other  pairs 
are  summarized  in  Table  3.  Overall,  significant  differences  are  not  found 
except  for  the  difference  between  group  1  and  group  4. 

The  resultant  estimates  from  (10)  can  also  be  applied  to  the  analyses  of 
ammunition  lots  which  are  not  used  in  the  experiment.  In  order  to  predict 
the  expected  quality  of  a  new  ammunition  lot  i  at  time  t  in  terms  of  depot 
and  vendor  characteristics  associated  with  lot  i  the  following  p)  is  used: 

fi  =  exp0o+$tt)/(l  +  exp0o  +$,$)  (12) 

where  fij  =  exp{\zt\+, ,  +7m2,'m)-  For  instance,  the  predicted  proportion 
defective  of  ammunition  lot  i  provided  by  vendor  A  after  5-year  storage  in 
depot  1  is 

p)  =  ea:p(-4.9034-(-ea;^(-0.9022)x5)/(l+Grp(-4.9034+ea;p(-0.9022)x5)) 

=0.0534.  (13) 

In  addition,  one  can  estimate  the  expected  time  t)  when  the  quality  of 
ammunition  lot  {  reaches  a  predetermined  level  p: 

%  =  [log{p/(l  -  p))  -  Po]/exp(yL*a+, .  +7m**n).  (14) 

Expected  time  t*,  when  the  quality  of  ammunition  lot  f  provided  by 
vendor  A  and  stored  in  depot  1  reaches  p=0.5  is  estimated  as 

t)  =  [log(0.5/(l  -  0.5))  +  4.9034]/ezp(-0.9022)  =  12.1.  (15) 
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For  the  four  groups  considered  in  this  example,  p}  when  fc=5  and  i*  for 
p  =  0.5  are  given  in  Table  4. 
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6.  SUMMARY  AND  DISCUSSION 

A  random  effect  logistic  regression  model  was  used  to  formulate  a  predic- 
tion model  for  ammunition  deterioration  rate  in  terms  of  depot  and  vendor 
information.  It  was  assumed  that  once  lots  are  selected  for  initial  inspec- 
tion, they  will  be  inspected  repeatedly  over  a  certain  time  period  without 
any  rectification.  The  proportion  defective  was  considered  a  measure  of 
ammunition  deterioration.  Regardless  of  different  sources  of  vendors,  the 
initial  quality  of  ammunition  was  assumed  to  be  homogeneous  as  was  its 
eventual  deterioration  with  the  passage  of  time.  The  log  odds  ratio  due  to 
unit  time  passage  was  defined  as  a  deterioration  rate  which  determines  the 
shape  of  the  logistic  curve.  Deterioration  rates  of  different  lots  of  ammuni- 
tion often  vary  and  their  variation  was  assumed  to  be  explained  partly  by 
a  function  of  depot  characteristics  and  vendor  information  associated  with 
each  lot. 

A  random  effect  logistic  regression  model  was  used  to  predict  the  dete- 
rioration rate  of  ammunition  in  terms  of  concurrent  characteristics  of  the 
ammunition  lot.  A  two-stage  analysis  was  employed  to  estimate  deterio- 
ration rates  and  relate  them  to  the  characteristics  involved  in  ammunition 
lot.  Using  the  approach  of  this  paper,  one  can  determine  the  appropri- 
ate time  for  reorder  or  renovation  of  ammunition  before  the  quality  readies 
unacceptable  levels  when  depot  and  vendor  information  associated  with  am- 
munition lot  is  available.  By  comparing  the  expected  deterioration  rates 
of  ammunition  stored  in  different  locations  one  can  improve  depot  mainte- 
nance policy  The  similar  approach  can  be  used  for  vendor  control. 

In  the  two-stage  model  used  in  this  paper,  we  considered  a  certain  type 
of  ammunition  and  its  deterioration  was  compared  in  terms  of  vendor  and 

18 


depot  characteristics.  In  order  to  compare  deterioration  rates  of  different 
types  of  ammunition,  the  following  approach  can  be  used:  First,  use  a  group 
of  relatively  homogeneous  types  of  ammunition  in  the  within-individual 
model,  and  secondly,  classify  them  not  only  in  terms  of  depot  condition 
and  vendor  information  but  also  by  types  of  ammunition. 

In  order  to  increase  the  accuracy  of  resultant  estimators  of  a  two-stage 
estimation,  sufficiently  large  number  of  inspection  which  captures  the  no- 
ticeable deterioration  of  ammunition  would  be  necessary  over  many  exper- 
imental lots.  However,  such  data  may  not  be  readily  available  in  reality. 
When  a  sufficient  number  of  repeated  measurements  of  actual  data  is  not 
available,  accelerated  testing  may  be  an  alternative  means  of  meeting  the 
time  and  budgetary  restrictions  of  the  experiment.  For  instance,  envi- 
ronmental tests  under  simulated  climatic  extremes  are  an  integral  part  of 
normal  engineering  tests.  These  tests  include  accelerated  depot  conditions 
at  the  test  site  when  the  characteristics  of  the  test  item  necessitate  the 
gathering  of  such  data.  Nonetheless,  when  only  a  small  number  of  data 
is  available,  resampling  procedure  such  as  bootstrap  method  (Dalai  et  al 
(1990))  can  be  used  to  make  inferences  concerning  deterioration  rates  at 
the  cost  of  heavy  computation. 

The  measure  of  rates  of  ammunition  deterioration  considered  in  this 
paper  was  the  proportion  of  defective  items.  When  other  variables  such  as 
muzzle  velocity  or  water  content  are  used  as  a  measure  of  deterioration, 
modifications  of  the  within-individual  model  may  be  necessary  to  take  into 
account  the  different  nature  of  the  measure.  However,  the  basic  concept  of 
a  two-stage  analysis  still  applies  to  any  modified  models. 

A  restriction  used  in  the  between-individual  model  was  that  a  lot  of 
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ammunition  would  be  stored  in  the  same  depot  in  the  experimental  period. 
A  prediction  model  that  reflects  changes  in  depot  conditions  associated  with 
a  lot  of  ammunition  is  under  investigation  and  is  left  for  further  research. 
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Figure  1.  Deteriorating  Patterns  of  Ammunition  Quality 


Table  1:   Number  of  Defectives  X^  Observed  in  Sample  Size  of  20  at  the 
j'th  Inspection  of  Lot  i 


group  1 


group  2 
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1 


lot  i 
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10 


1 
2 
3 
4 
5 
6 
7 
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0 
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na 


0 
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1 
1 
1 
2 
3 
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Pi     0.3853     0.4490       0.4172     0.3671     0.3853 
r~     0.0071     0.0055       0.0063     0.0077     0.0071 


0.4446       0.3649     0.4969 
0.0033       0.0045     0.0028 


0.4446     0.2310 
0.0033     0.0088 
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10 
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6 
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5 

9 

7 

5 

6 

5 

13 

10 

8 

8 

10 

7 

10 

na 

na 

na 

na 

na 

13 

10 

10 

13 

9 

Pi 

0.4884 

0.4355 

0.4535 

0.4143 

0.6290 

0.5582 

0.5009 

0.5009 

0.5536 

0.4776 

0.0019 

0.0022 

0.0021 

0.0023 

0.0018 

0.0014 

0.0014 

0.0014 

0.0014 

0.0014 

/ 


Table  2:  Fitting  the  Between-Individual  Model 
ji     standard  error  (^) 


1  -0.9022 

2  -0.8336 

3  -0.7184 

4  -0.6567 


0.084186 
0.068770 
0.052690 
0.047560 


Table  3:  95%  Confidence  Intervals  for  the  Difference  between  Expected  Log 

Deterioration  Rates  of  Two  Groups 

group  /  (vendor  and  depot)     -     group  /'  (vendor  and  depot)     (lower  limit,  upper  limit) 


1 

(A  and  1) 

1 

(A  and  1) 

1 

(A  and  1) 

2 

(B  and  1) 

2 

(A  and  1) 

3 

(A  and  2) 

2  (B  and  1 

3  (A  and  2 

4  (B  and  2 

3  (A  and  2 

4  (B  and  2 
4  (B  and  2 


(-0.2990,  0.1618) 
(-0.3943,  0.0267) 
(-0.4504,  -0.0406) 
(-0.2989,  0.0685) 
(-0.3542,  0.0004) 
(-0.2122,  0.0888) 


Table  4:  Predicted  Proportion  Defective  at  year  5  and  Estimated  Time  for 
p  =  0.5 


» (%) 

t)   (year) 

5.34 

12.1 

6.12 

11.3 

7.83 

10.1 

9.02 

9.5 
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